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Abstract 

We obtain sufficient conditions for the stability of the synchronized 
solution for certain classes of coupled dynamical systems. This is ac- 
complished by finding an analytic expression for the transverse Lya- 
punov exponent through spectral analysis. We then indicate some 
applications to population dynamics. 

1 Introduction 

The study of coupled dynamical systems has received considerable atten- 
tion recently for its interest from the mathematical, physical and biological 
point of view, see for instance ^F3\, 0, [U, [IS], jUJ and among oth- 
ers. One concern about such systems is to whether or not they will present 
synchronization phenomena and to whether or not such synchronization is 
stable. 

The system below describes the evolution of a system consisting of d 
identical subsystems where, on every iteration, each subsystem undergoes 
its common local evolution determined by /, followed by a density dependent 
coupling process encoded by C and tp. The system can model a population 
consisting of d patches, Xj, j = 1, . . . , d, where, in the absence of migration, 
patch j is controlled by a local dynamics x\ +l = f{x° t ). When migration is 
present, ip(f(x J t )) individuals leave patch j and are distributed with density 
Cji on patch i. The global dynamics is then 



j" 



if i = f(4) - vtfte)) + Yl c Mf(4)y, i,j = i,"',d. 
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Here / is a C 1 map on [0, oo), C = [cij] is doubly stochastic, that is, c%j > 

d d 

and y ] Cjj = ^Jcy = 1- Furthermore we will assume ip differentiable 

i=l j=l 

a. e. with bounded. The models in population dynamics considered in 
G] EH [H] H^l are au particular cases of © for special choices of 

C and ip. 

The condition on C being doubly stochastic reflects that there are no 
losses during the migration process. It is also necessary for the invariance 
of the diagonal of the phase space, that is, for x\ = x\ = xt to be a solution 
of (pQ), where each x\ satisfies x J t+1 = f(x J t ). 

In this paper we obtain sufficient conditions for the stability of the afore- 
mentioned synchronized solutions. The criteria involve the Lyapunov expo- 
nents of the one dimensional map / and of the codimension one transverse 
dynamical system to be defined below. 

The paper is organized as follows. In the next section we provide a 
criterion for stability for general systems. In §3 we improve our result for 
the case of normal operators. In §4, we consider a system where the cou- 
pling/migration process is time dependent. We formulate and prove the 
corresponding results of the previous sections in this setting. In the last 
section we indicate some applications to population dynamics. 

Previous results treat only the case where <p' is a constant, or a 2-valued 
step function as well as particular examples of matrices C and in these cases 
we recover such results. Our treatment only requires <p' to be bounded and 
C to be doubly stochastic and irreducible. This last condition is easily seen 
to be necessary, see below. Moreover we are unaware of any treatment of 
the systems in §4 in the literature. Thus we extend some of the results of 
|. 0; HH1> and others to these more general situations. 



2 Stability: General case 

In order to understand the behavior of orbits starting at nearby points of 
the diagonal of the phase space, we first linearize If Jt = [ajj] denotes 
the Jacobian matrix of (^Q) restricted to the synchronized orbit, we have 

f'ixt) (1 - (1 - cuW (f(x t ))) , for i=j 
f\x t )ip'{f{xt))cij, for i^j 

We have J t = f'(x t )H t , where H t = I — <p'(f(x t ))B, where B = I-C. 

We will assume that C is irreducible. This is almost a necessary con- 
dition for otherwise it permits the existence of uncoupled unsynchronized 
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subsystems that are each synchronized. In this case we can apply Frdbenius 
theorem |1U| to show that A = 1 is the simple dominant eigenvalue of C, as- 
sociated to the eigenvector u = (l,...,l). This furnishes the decomposition 
M d = W, where W is C— invariant d — 1 dimensional subspace. Under 
these conditions 

" 

A 



B = P- 1 



P 



where P is the matrix of the appropriate change of basis. This decomposition 
implies that the stability of the synchronized solution of Q is a consequence 
of the stability of the trivial solution of the transversal component, wt, which 
satisfies 

w t+1 = f'(x t ){l-<p'(f(x t ))A)w t . (3) 

We will show that under a certain integrability condition the map above 
is in fact a contraction which in turn implies the stability of wt = 0. The 
analysis of (J3J) will be based on the Lyapunov exponents (see ^H], ^H]) of 
©. Define 

n-l 

K n {x) = [J f{f k {x)) (l - v'U k+1 {x))A) 

k=0 

where f°(x) = x and f k {x) = f(f h ^ 1 {x)) for k > 0. Clearly if fC = 
limsup || K n || l / n satisfies /C < 1 we have that © is a contraction. Now 
observe 

K n {x) = L n {x)K n {x)= m /'(A*)) J (l[l-ip'(f k+1 (x))A\ (4) 

\fc=0 / \k=0 / 

L n {x) depends only on the local dynamics / while A n (x) reflects also the 
effects of ip and C. Let p be an invariant measure of the local system. 
Define for x > 0, ln + (x) = max(ln(x), 0). By Birkhoff's ergodic theo- 
rem, if ln + \f'\ G L x (p), there exists lim n ^ J2k=o ^ n \f'(f h ( x ))\ f° r /°" a - e - 
x. For p ergodic, this limit, call it L, is independent of x and is given by 
J °° In |/'(s)| dp(s). L is the Lyapunov exponent of the local system governed 
by/- ' 

Similarly the ergodic theorem of Oseledec |15j implies that if 



/•oo 

/ \n + \\I-^(s)A\ \ dp(s) < oo, 

J 



there exists lim n — In ||A n (a;)|| =: lnA(x) for p-a.e. x and this limit is inde- 
pendent of x provided p is ergodic. Our first result is 
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Theorem 1 Consider the system ^) where f is a C 1 map, p> bounded, C 
doubly stochastic, irreducible such that\n + \\I —Lp' (s)A\\ andln + are in 

L l (p), where p is an f —invariant measure. Let L = sup,,, lim n \L n (x)\ 1 ' n and 
A = sup x lim n ||A n (x)|| 1 / n . If L A < 1, there exists a set E with p(E) = 1 
such that for all x G E, the synchronized solution of 0j is asymptotically 
stable. 

Proof: By Oseledec's theorem there exists a set E with p(E) = 1 such that 
for all x in E, lim n ||A n (x)|| 1//n = A(x). We claim that for all x in E 



A{x) < exp (J ln+ 1 1 1 - <p'(s)A\\ dp(s)) 



(5) 



By the continuity of the norm and the function ln(-), the dominated conver- 
gence theorem implies that we only need to prove (J5J) for ip' simple. 

Let ip'(x) = Yjk=l a kXE k (x), where E k are measurable and disjoint with 
p(E k ) > and such E c[j k E k . For 1 < k < r, let 

»{0 < j < n : P(x) E E k } 

Pk,n = • 

n 

This gives 

A n (x) = Y[(! ~ a k A) npk - n , which imply 
k 

||A n (x)|| < H\\I-a k A\\ n P^, and 

k 

ln\\A n (x)\\ < £ 

n Pk,n 

ln||/-a fc ,4|| 

k 

Birkhoff 's ergodic theorem applied to \E k shows that for 1 < k < r, we have 

limp kn = p(E k ). 

n 

Therefore, given e > 0, there exists no such that for n > no, we have 
Pk,n < (1 + e)p(E k ). Since ln(-) < ln+(-), 

-i r roo 

-ln||A n (x)|| < (1+e) Vln+ \\I-a k A\\ p(E k ) = (1+e) / ln+ \\I-ip'{s)A\\ dp{s) 

Note that the right hand side is independent of x. Thus we have 
||^n(*)H 1/n < \L n (x)\ 1/n \\K(x)\\^ n < (l + e)LA 1+s 
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and since e > is arbitrary we have 

lim||K n (x)|| 1/n < L A < 1 

n 

and thus the transversal map © is a contraction. □ 
2.1 Remarks: 

(i) if E = [0, oo), L A < 1 is necessary for the stability of the synchronized 

solution. 

(ii) in certain cases, which include the case A semi simple, the norm 

— ip'(f k (x))A\\ is simply the spectral radius, <j-i(H t ), of the restric- 
tion of B to the subspace W . Therefore, with the same hypothesis on 
/, ip and C we have 

r-oo 

A<A 1= / cr_i(i?^ (s) ) dp(s) 
J o 

and thus LAi < 1 is a sufficient condition for the stability of synchro- 
nized solution of (^Q). 



3 Stability: normal operators 

In the case of normal operators, that is, AA* = A* A, one can improve 
the previous result with the help of the functional calculus of ^1] and [2D] . 
The proof of Oseledec's theorem in shows that lim n (A* A n (x)) ^ = 
A(x) in operator norm. Spectral analysis of A(x) determines the Lyapunov 
exponents and respective subspaces. Even though in our case C represents a 
dx d matrix, the discussion below is valid for any bounded normal operator 
in a Hilbert space. We present the results in this generality since it will be 
applied in the future to systems more general than We will show that 
for A normal, we have 

A(x) = exp ^ In \I - tp'(s)A\ dp(s)j (6) 

Without loss of generality we can assume that for s in E, I — ip'(s)A is 
nonsingular. This implies that for all A G &(A), In |1 — \(p'(s)\ is continuous 
on E. The spectral theorem for bounded normal operators (Theorem 12.23 
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in [10]), allow us to define \I — tp'(s)A\ and In \I — \(p'(s)A\ by 

A = [ XdP(X) 

Ja(A) 

In \I - <p'(s)A\ = [ In |1 - X<p'(s)\ dP(X) 
Ja(A) 

where dP(X) are the spectral projections associated with A. Our hypothesis 
on the spectrum of / — (p'{s)A together with Fubini's theorem imply 

/ \n\I-tp\s)A\dp(s)= / / ln|l- \(p'(a)\ dp(s) ) dP{X) 

JO Ja(A) \Jo J 

and since this last integral defines a bounded operator, we take exponentials 
to obtain 



exp 



jf In \I - <p'(s)A\ dp(s)^ = exp \J ^ In |1 - Xtp'(s)\ dp(s)^j cLP(A)^ 

and thus the right hand side of © is well defined. The spectral mapping 
theorem then implies 

cr^expQ! ln\L-ip'(s)A\ dp(s))j = jexp (J In |1 - Xip'(s)\ dp(s)^j : X € o{A) 
Our next result is 

Theorem 2 Let f, ip 1 , C, and L be as in Theorem^ Let A be the spec- 
tral radius of Then, if L A < 1, the synchronized solution of 0j is 
asymptotically stable. 

Proof: The above considerations imply that in order to prove the theorem 
we only need to establish ©. We will do so first assuming <p' simple. Then a 
standard limit process extends the result to ip' in L°°. Let ip' = Y2k=i PkXEk 
where E k are disjoint measurable. Note that since p' is bounded, the integral 

m loD is well denned. As before we define Pkn = • 

w n 
Then 

N / N \ 

(A;A n )^ = Y[\L -<p k A\P k > n = exp ( ^2 Pk,n In \I — VkA\ J 

fc=i Vfc=i / 
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Since for all k, lim n p nj fc = p{Ek), the continuity of exp(.) and the Lebesgue 
dominated convergence theorem imply 

lim exp (^2 Pk ' n ln 1^ ~ VkA^j = exp ^ J In 1 1 - tp'(s)A\ dp(s)^ . 

proving ©. This finishes the proof of Theorem |2JD 
3.1 Remarks: 

(i) Since for all A G a (A), |1 - \(f'(s)\ < \\I - (p'(s)A\\, we have A < 

exp f J ln + 1 1 J — ip'(s)A\ \ dp(s)^j , therefore Theorem|2]is an improve- 
ment of Theorem ^ 

(ii) Note that all the effects of the weighted network relevant to synchro- 

nization are reflected on the spectrum of A and can therefore, in certain 
cases, be independent of d. It would be interesting to extend such re- 
sults to the d = oo case since Frobenius theorem as well as Oseledec's 
ergodic theorem are not valid without further assumptions. 



4 Extensions 



The analysis above can be applied to more general systems. One such in- 
stance is the case below where the coupling/migration process no longer 
depends on the density but instead obeys a seasonal dynamics, that is, on 
each cycle a time dependent fraction pt of each patch will mix according to 
a time dependent distribution Ct = [c*-] as follows 

d 

x 3 t+1 = f{x{)- p t fixD + ^c^pt f{x\)] i,j = l,---,d. (7) 

i=i 

Again / is a C 1 map on [0, oo). For g continuous on [0, 1] and po arbitrary, we 
assume that p t = p t (po) is given by p t+ i = g{p t ) ■ Similarly let {C s } s6 [ 0j i] 
be a family of doubly stochastic irreducible matrices. If h is a continuous 
map on [0,1] and so is arbitrary, define st+i = h(st) and Ct = Ct(so) by 
Ct+i = Ch( St y Assume that G(p, s) = p C s is a measurable operator valued 
map with respect to the product measure v x rj on [0,1] x [0,1], where 
v and r] are, respectively, a ^-invariant and /i-invariant ergodic measures 
on [0,1]. Under these assumptions the diagonal of the phase space, x\ = 
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x t+i = x ti 1S a synchronized solution and we are interested in its stability. 
This system falls under the general theory of Random Dynamical Systems as 
presented in [Q and thus obey, under the appropriate integrability condition, 
a mutiplicative ergodic theorem. Our approach is direct and avoid the use 
of this theory. Moreover, due to the particular nature of (0), we are able, 
as in the previous sections, to obtain much more precise information about 
the limit operators than is given by the ergodic theorem alone. 
The Jacobian matrix of (J7J), Jt = \oi\A is now given by 



yielding J t = f'(xt)(I — fitBt) where B t = I — Ct- Our first task is to 
decompose the linearized system into diagonal and transversal components 
and then estimate the respective Lyapunov exponents by the appropriate 
ergodic theorem. This cannot be done in general for arbitrary families {Ct}, 
therefore we will consider special cases of increasing generality. 

4.1 Simultaneous Diagonalizable Matrices 

This include the families of commuting symmetric matrices and of circulant 
matrices. In this case there exists a matrix P such that for all t, Bt = 
P^ 1 Mt P where Mt is a diagonal matrix with entries {1, A2(i), • • • , Xd(t)} in 
its diagonal. Each Xj(t) is measurable and satisfies |Aj(t)| < 1 and Xj(t) ^ 1 
by Frobenius theorem. In the symmetric (commuting) case the A,-(f) are real 
and lie in (—1, 1). We obtain that the transversal component of (jZJ) satisfies 



where Dt is the (d — 1) dimensional {Xz(t), . . . , Xd(t)} diagonal matrix. As 
before, the synchronized solution of ® is stable if and only if wt = is a 
stable solution of (JHJ). For each x, fio, so in [0, 1], define L n {x) as in Q), and 
for j = 2, . . . , d, define 




wt+i = f'(x t ) (I - fit D t ) w t . 



(8) 



n-l 



Ajn(W), so) = ]J I 1 -/(^o)A i (/i fc (s ))| 



fe=0 



The analogue of Theorem |^1 is 



Theorem 3 Consider the system Q) under the above conditions. Assume 
in addition 



S 



(i) ln + |/'(s)| belongs to L l (p), where p is a f— invariant measure on [0, oo). 

(ii) For all 2 < j < d, ln + [1 — p\j(s)\ belong to L l {y x rf). 

Then there exist sets E C [0, oo) and F C [0, 1] x [0, 1] with p(E) = [y x 
rf) (F) = 1, such that for all x G E and (po,so) G F, the limits L{x) = 
lim n \L n (x)\ 1 ^ n and Aj(no, sq) = lim ra \Aj n (/Uo, sq)^^ exist. Moreover if L = 
sup x L(x) and A = sup JjA10 S0 Aj(po, so) satisfy L A < 1, the synchronized 
solution of ^ is asymptotically stable. 

Proof: The existence of L(x) for x in a set of full measure follows, as before, 
from Birkhoff's ergodic theorem applied to lnL n (-) and (i). In addition 
the limit is independent of x provided p is ergodic. The decomposition Q 
allow us to derive the transversal component directly without making use of 
Oseledec's theorem. For 2 < j < d, (po, sq) G [0, 1] x [0, 1], write 

A in (> ,so) = exp ( - 5^1n|l - g k (p ) ^j(h k (s ))\ J 
\ n fc=0 / 

Condition (ii) and a (d — l)-fold application of Birkhoff's ergodic theorem 
imply the existence of a set F, with v x rj (F) = 1 such that for all j, all 
(po, sq) G F, there exists Aj(po, sq) = lim n Aj n (po, sq). If, in addition, v x r/ 
is ergodic, the limit is independent of (po, sq) and is given by 

Aj = exp / In 1 1 - p\j(s)\ d(y x rj)(p,s) (9) 

This gives diag{A2, ■ ■ ■ , A^} as the analogue of ©• Its spectral radius is then 
A = maxj Aj. Thus, if LA < 1, the synchronized solution is asymptotically 
stable. □ 



4.2 Symmetric (non commuting) Matrices 

In this case we no longer have a decomposition yielding a diagonal transver- 
sal component, like (jSJ), and therefore our conclusions are somewhat weaker 
than the previous theorem. Since each Ct is symmetric, doubly stochas- 
tic and irreducible, A = 1 is the dominant eigenvalue with corresponding 
eigenvector v = {1, . . . , 1}. Moreover the {d — 1)— dimensional subspace 
W = (Kv) 1 - is invariant for all Ct- This provides a decomposition like © 



B t = P 



-i 







At 



P 



(10) 
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Accordingly the transversal component of (J7J) now satisfies 

w t+1 = f'(xt) (I - lit A t ) w t (11) 
The following theorem is the analogue of Theorem ^ in the present situation 



Theorem 4 Consider the system Q) where f and fa are as * n Theorem^ 
and assume {Ct} as above. For each (/io> s o)> define 

n-l 



k=0 



Assume that ln + \\I — fiA h ^\\ belongs to L l {y x rj). Then there exists a set 
F C [0, 1] x [0, 1] of full v x r]-measure such that for all (fj,Q, so) 6 F , the limit 
A(/i ,s ) = lim n ||A n (^ ,s )|| 1/n exists. Moreover if A = sup (w))S0 ) A(/i ,s ) 
satisfies LA < 1, the synchronized solution of Q) is asymptotically stable. 

Proof: Consider the map Q on [0, 1] x [0, 1], given by Q(fJ., s) = (g(fJ-), h(s)). 
Clearly v x r\ is Q invariant and by assumption, A n (ji,s) above defines a 
(v x 7/)-measurable cocycle. The existence of A(/i ,so) f° r (/^o^o) on a 
set of full measure then follows from Oseledec's ergodic theorem applied to 
A n (/i,s). We also have that for a.e. (fio,so) 

A(// , s ) < exp ( / ln+ \\I - fiA h(s) \\ d(y x r/)(/x,s) (12) 

\J[0,l]x[0,l] / 

The proof of ()12j) follows the same path as the proof of (J5J) , thus we will omit 
the details. Note that since each At is symmetric, there exist \{t) £ o~(At) 
for which \\I — n A t \\ = \l — n\(t)\. X(t) is the lowest eigenvalue of Ct- If, 
in addition, the measure v x n is ergodic then A(fiQ,so) is independent of 
(Ho, sq) and in this case we have 



A = exp / In |1 — /it X(h(s))\ d(v X r))(fi,s)\ 

V-/[0,l]x[0,l] ' ' / 

as a consequence of Birkhoff's theorem. As in the previous theorems, the 
condition L A < 1 implies the asymptotic stability of the synchronized solu- 
tion. □ 
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4.3 Remarks: 

(i) The assumptions on the continuity of g and h above can be relaxed. 

All that is needed is that each map possesses an invariant probability 
measure, v and rj. 

(ii) Similarly to the density dependent migration, the results above do not 

extend to the d = oo case even though once again the effect of d is 
encoded in the joint spectrum of {Cj}. 



4.4 Examples: 

Some special cases of © and (|12|) are worth mentioning. 

(i) if so is a periodic point for h, that is, h p (so) = sq, and we take n = 

I Efc= ^ fc («o)' then © becomes 

Aj=exp\-y~] [ ln|l -fxXj(s k )\ du(fi) ] 
V t^o J M J 

A similar expression corresponding to the case of a periodic point for 
g also holds. 

(ii) if g = h, and rj = i>, but [Aq ^ so, we get for @ 

Aj = exp / In 1 1 — [i Xj (s) | d(v X v) (/x, s) 

V-A ' 1 !^ ' 1 ] ' J 

(iii) if in (ii) we have, in addition, /xo = sq, that is, the term \i t A t in (fTT|) 
is of the form g t (so) A g tr SQ -\, then Q becomes 



Aj = exp / In 1 1 — s Xj(s)\ dv{s) 



, '[0,1] 

Similar corresponding expressions can be obtained, without difficulty, to 
represent (|12|) in the special cases above. We leave the details as well as the 
formulation of other special instances of (|12|) to the interested reader. These 
formulas are useful in cases where v and r\ are known, say Lebesgue, for they 
permit the exact calculations of the Lyapunov numbers of the transversal 
dynamical system. 
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5 Applications 



In [3j we apply the above results to several special cases of interest in pop- 
ulation dynamics. In this section we will restrict ourselves to the density 
dependent system (^Q) and leave the corresponding formulations of the results 
related to the time dependent system (JJJ) to the interested reader. 

A direct problem consists of determining stability once /, if and C are 
given. In such situations we compute the eigenvalues Ai, . . . , A^ of C. Dis- 
carding the eigenvalue 1, we evaluate the integrals defining the spectrum of 
A(x) above by Birkhoff 's ergodic theorem 



then the maximum of these values gives the spectral radius A and we can 
then determine stability once the Lyapunov exponent for / is known. We 
note that in order to use (|13|) we need the /—invariant measure, p, to be 
absolutely continuous with respect to Lebesgue measure. Such measures, 
called physical or SRB (Sinai-Ruelle-Bowen) measures, are known to exist 
in certain important cases, such as expansive maps, piecewise monotonic 
maps of the interval as well as Axiom A diffeomorphisms, even though a 
recent result of Bochi-Yoccoz shows that they are not typical in the C 1 - 
topology. Further information can be found in |2(ij . |17j . |18j and |15j . 

One feature of Theorems 1 and 2 is their continuous dependence on <p' ', 
that is, C 1 -uniform topology. This follows directly from the formulas for 
A above. In [3] we give examples that show that this dependence is not 
continuous in ip in the C°-topology. 

We can also make use of the above results to study an inverse problem. 
In this case / and <p are given and one is interested in finding a double 
stochastic matrix C, if possible, yielding the desired stability behavior for 
the synchronized solution of (JJ). In this regard, Theorem H3 below gives a 
partial result for symmetric matrices. We will make use of the following 
result from |12j . 

Proposition 5 The following are corollaries 7 and 8 from llty . 

(i) If {A2, • • • , Xd} 6 [—l/(d— 1),1], then there exists a symmetric doubly 

stochastic matrix C with <r(C) = {1, A2, • • • , A^}. 

(ii) If X € (—1, 1], there exists a positive symmetric doubly stochastic matrix 

C such that X £ cr(C). 




k=0 



(13) 
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The proof of Proposition |S] in ^2] provides algorithms to find the matrices 
C in (i) and (ii) above. 



Theorem 6 Let f and tp be as in Theorems and Let L denote the 
Lyapunov exponent of the one dimensional map f. For A E [0,2], define 



F(X) = exp / In 1 1 - \<p'(s)\ dp(s) , m = inf F(X) and M = supF(A). 



(i) if Lin > 1 the synchronized solution of 0) is unstable for all symmetric 

configurations C . 

(ii) if L M < 1, the synchronized solution of 0j is stable for all symmetric 

configurations C . 

(iii) if L G (l7>^); ^ * s possible to find a symmetric doubly stochastic 
matrix C such that the synchronized solution of 0j has a prescribed 
stability behavior. 

Proof: Parts (i) and (ii) are direct consequences of Theorem [2 and the 
fact that F(\) gives the spectrum of F{A). To prove (iii) we first note 
that any non negative symmetric doubly stochastic matrix has its spectrum 
contained in (—1, 1]. This follows from Gershgorin's theorem [Jj which states 
that cx(C) is contained in the set {A G C : Vi, |A — ca\ < \ c ijW- The 

hypothesis on C then imply that o~(C) C [—1, 1]. If C is positive, it is not 
difficult to show that —1 cannot be an eigenvalue of C. Since B = L — C , 
we have cr(B) C [0,2). Theorem H3 now follows easily from Theorem |^1 and 
Proposition E|D 
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